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1. INTRODUCTION 

In the last few decades, the dangers of fossil fuel exhaustion have led research centers and 
companies involved in the field of renewable energy to offer more affordable and less expensive technical 
solutions. Among these solutions, the doubly fed induction generators (DFIG) which is become extensively 
employed in the wind power generation systems. This is due to their remarkable characteristics compared to 
the permanent magnet synchronous generator (PMSG) [1]. Some of these characteristics are the adaptability 
with variable wind speed by simply changing the rotor operating frequency and the regulation of the active 
and reactive power injected into the grid. Also, the converter related to the rotor is designed to support energy 
exchange rates between 25% and 30% of the generator's global energy [2]. Many research works have 
addressed the control of DFIG, the most known strategies are the field-oriented control (FOC), the direct 
torque control (DTC) and the direct power control (DPC) [3], [4]. 

The first method has some difficulties, among them, the necessity to use coordinate transformations, 
the current PI regulators must be tuned and they depend strongly on DFIG parameters, in other words, this 
method is not robust. The DTC strategy has been introduced in the aim to overcome these problems [5]. It is 
known for being performing, robust and less sensitive to parameters variations of the DFIG. Unfortunately, 
the high ripple in developed torque and power, in addition to the non-constant switching frequency led to a 
non-linear behavior in the torque and the power hysteresis comparators. The progress of power electronics 
permitted the emergence of rapid switches, this has enabled the development of new control strategies such 
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as the direct power control (DPC) Figure 1 [12], [3]. This method has mitigated considerably the problems 
mentioned. 





Figure 1. Schematic diagram of the DFIG-based wind energy generation system 


The existing DPC strategies are classified into two major categories; the voltage-based DPC and 
virtual flux-based DPC [6]. Figure 2 gives more details about the different sub-classes of these control 
strategies. Our proposal aims to analyze and enhance the performances of classical DPC. Also, the 
assessment of the DFIG response to a variable wind profile [7]. 
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Figure 2. Block diagram of DPC strategies 


This paper aims generally the study of the direct power control strategy (reactive-active) applied to 
the rotor-side converter (RSC). The main contributions of this work are, a) design of direct power control 
strategy that ensures a high level of performances, b) guarantee the extraction of the maximum power and 
reaching a unit power factor by a decoupled control of active and reactive powers, and c) applying a fixed 
and aleatory wind profile to wind energy conversion system to assess its capability in different usage 
situations. 

The remainder of this article is organized as follows: First, Section 2 presents the model of the 


Performance analysis and enhancement of direct power control of DFIG based ... (Mohamed Amine Beniss) 


1036 O ISSN: 2088-8694 


studied wind energy conversion system. Then, section 3 describes the adopted DPC strategy. After that, 
Section 4 discusses the simulation results. Finally, section 5 provides some conclusions of this work. 


2. MODEL OF THE WIND ENERGY CONVERSION SYSTEM 
2.1. Mathematical model of the wind turbine 
The wind turbine aerodynamic input power can be expressed as [8]-[10]. 


P, => pSV (1) 


where p refer to the air density, while S represents the area of the wind wheel, V denotes wind speed. 


The formula expresses the mechanical power gained by the turbine from the aerodynamic input 
power. 


Pe = Cp As BY By = 5" CAB) pSV" (2) 


Where the coefficient Cp represents the efficiency of the wind turbine blades. It is a function of À and B which 
refer respectively to the tip speed ratio and the blade pitch angle. The parameter à can be obtained by 
using (3). 





V (3) 


With: 
R: the blade radius. 
O: The angular speed of the turbine. 


The function Cp can be approximated by this well-known expression [1], [8]. 


S Caine a LE (4) 
C, =0.5-0.0167 -(B—2)) sin 55-034 B- 


The coefficient Cp reaches it maximum value Cpmax=0.4998 when B=0 degree and A=8.095. Based 
on this point one can define the rotor speed that corresponds to the (MPPT) [9]. Figure 3 depicts the time 
variation of the power coefficient Cp. This result has been obtained after the simulation of our system, the 
simulation curve of Cp Figure 3 proves the robustness of our MPPT algorithm. The Cp maximum value is 
maintained at about 0.4985 despite the variation of the wind speed [8], [11]. 


Power coefficient 





Time(s) 


Figure 3. Power coefficient (C, ) 


2.2. DFIG modeling 
Many research works have been devoted to the modeling of wind power generators. Especially the 
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DFIG, in order to simplify the control of the DFIG we have to model the system in the two-phase reference 
Figure 4 [12]. By using Park transformation, we convert the three-phase representation of the system to 
another one in the d-q frame. Thus, the overall electrical model can be obtained. The explain this clearly [3], 
[13], [14]. 

— Voltages across the stator terminals 


. _ dy, 
Va = R, la + Ge ~ Qo, W sq (5) 





dV, 
Va = R, ds, + T + Q, Wa 


— Voltages across the rotor terminals 





V,=R.i,, + ~O,Wy (6) 
Va =R, it Ku +O.YW,, 
With: 0.-@.=@,, 
The basic magnetic equations related to the flux expressions in the reference frame (d, q) are given 
as. 
— Flux at the stator 


W sa = Liga gä M i, (7) 
W sa E L, Lg + M i; 


— Flux at the rotor 


Wa = Leha + Mda (9) 
Y,, =Li1,,+M i, 


The electromagnetic torque expression as a function of the rotor flux and the rotor currents is given 
by [15]. 


3 
Con = 5 Pla Wiad — La Y) (9) 
By using the fundamental equation of dynamics, we deduce [16]. 
E E ESS E (10) 
dt 


Where Cm, fy and J are respectively the mechanical torque applied to the DFIG, the viscous friction and the 
wind turbine inertia. 





Figure 4. d-q equivalent circuit of the DFIG in synchronous coordinates 
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2.3. Voltage converter modeling 

In this subsection, the IGBT transistors are considered as ideal switches. Consequently, one can 
derive the relationship between control configurations and the electrical voltages of the AC and DC sides of 
the inverter. The general topology of the two-level voltage source inverter (VSI) is shown in Figure 5 each 
inverter arm contains two electronic switches IGBT [25]. 

All the information mentioned are summarized in the mathematical model [17]-[19]. 


y [2 = “HS 
==S|-1 2 -1]/5, 
A et as 


Where, Var, Vor, and Ver are are the output voltages of the inverter. 
It is underlined here that the switches states can be represented by three Boolean variables such as, 
a) Si = 1 if T; is closed and T;is open, and b) Si = 0 if T; is open and T; is closed. With (i = 1,2,3). 





Vy 
Vy 
ya 








Figure 5. Simplified topology of the two-level VSI 


3. DPC CONTROL STRATEGY 
3.1. DPC principal 

The DPC strategy aims to directly control the stator active and reactive power of the DFIG [24]. It is 
based on the change of the rotor flux vector in the d-q coordinates. This can be done by controlling the 
voltage vector at the output of the rotor side converter (RSC). As mentioned before the control of IGBTs 
generates the desired voltage vector [13], [20]. 

To select the appropriate state of the converter's switches. The instantaneous sector zone of the 
rotor’s vector flux, the error between the real active power and the reference one and the error between the 
estimated reactive power and its reference must be calculated. Two hysteresis comparators are used to 
guarantee that both powers are confined between the predefined upper and the lower bounds. These operating 
steps are illustrated more clearly in Figure 6 [21]. 


Estimation block 





Figure 6. Detailed scheme of the DPC strategy 
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3.2. Active and reactive powers formulation 
The expression of the rotor vector voltage can be formulated as [13], [17]. 


V. = {/2U.(S,+S, e" +8, e") (11) 


Figure 7 displays all possible orientations of the voltage applied to the rotor. It is noticed that there are eight 
different voltage possibilities two of them are zero vectors (Vp and V3). 

At this stage, we establish the mathematical relationship between rotor flux and voltage in the d-q 
reference frame. 


FoR 4 (12) 
dt 





For the sake of simplicity, we neglect the rotor resistance for the high and medium power machines (6). 


. day 
Ves (13) 
” dt 


Knowing that both simulation platforms and electronic calculation units work in discrete time. The (6) needs 
to be discretized. Therefore, applying this process to (6) yields. 


i 
[V =y; fin- 7 Gnt) (14) 
0 

We obtain: 
Wi (fin) (int) = V7, (15) 


The (16) and (17) give the expression of stator active and reactive powers in the d-q reference frame, 
respectively. 


P =5 (Vi Ia FV pdt ca) (16) 
= Gh E, (17) 
= 56 sd *~sd "sd a 


After the development and simplification of the equations mentioned before, we can obtain the 
results [13]. 
P, =K, |y, |sin (ô) (18) 


Q. =K, | K, —|y,|cos(6) | (19) 





y,|sin(d) 
and lW7,|cos(d ) will influence P andQ.. Figure 8 graphically explains the equations. It clearly shows how 





K,, K,and K, are constant parameters of appropriate dimensions. It is remarkable that changing 


the rotor voltage vector can modify the direction and magnitude of the rotor flux. The module and the phase 
of the rotor flux in the d-q reference frame are provided by. 


0, = arean 2 (19) 
Wig 


Y, = [Vu Wa (20) 
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Figure 7. Voltage vectors delivered by the Figure 8. Vector diagram describing some possibilities 
inverter for changing the location of the rotor flux 


3.3. Hysteresis power regulators 

Many solutions offer the possibility of maintaining very efficient operation of the DFIG, but among 
them, the use of two or three-level hysteresis regulators remains the most appropriate choice for controlling 
the active and reactive power levels. By inserting these comparators in the DPC diagram, we guarantee that 
the exchange of active and reactive power between the DFIG and the grid follows the reference values. 
Hence, the selection of a rotor voltage vector must reduce errors between the measured and the reference 
powers (active and reactive) and confine them between the hysteresis upper and lower bounds [22] [9]. To 
maintain Q, between the desirable bounds, we choose a two-level hysteresis regulator as depicted in Figure 9. 
The Boolean function So increase Qs when its value is equal “1” and decrease Qs otherwise. Here are the 
instructions that control the behavior of the hysteresis regulator at two levels. 


if AQ > £o > Sọ =1 
if -£o < AQ < sand = > 0 = S, =0 (21) 


A 
if -£o £ AQ < egand “= <0 = Sy =1 


if AQ < —£, => So =0 


While to control the active power Ps we use Three-level hysteresis regulator as displayed by Figure 10. In 
this case, the function Sp can take three values: if “1” increase P ,”-1” decrease P., and “0” to hold Ps. The 


algorithm contains the instructions that govern the functioning of the three-level: 


if AP =E, >S, =1 
AP 
ifO< AP < Ep and ae eee 
t 


if O < AP < €, and Ag <0=5,=1 
t 


(22) 


if AP ~—€, => Sp = -1 
AP 
f =E, <AP<O and Pa >O> S, =0 
t 


AP 
if —E, S AP <O and a nds aa 


AQ 





Figure 9. Hysteresis reactive power controller. Figure 10. Hysteresis active power controller. 
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3.4. DPC Switching table 
Based on instantaneous input values such as the rotor flux vector sector and the state of the 
hysteresis output signals the switching table selects the next rotor voltage vector to get the estimated rotor 


flux (Table 1). Each rotor voltage vector corresponds to a specified switches state combination of the RSC 
[23]. 


Table 1. Switching table for DPC method of DFIG 


Sector Oy 
1 2 3 4 5 6 
[330,30] [80,90] [90,150] [150,210] [210,270] [270,330] 
Qs Ps Voltage vector 
1 V; Ve Vi V2 V3 Va 
1 0 V3 V3 V3 V3 V3 V3 
-1 V3 V, Vs V6 Vi V> 
1 Ve Vi V2 V3 Va V; 
0 0 V3 V3 V3 V3 V3 V3 
-1 V2 V3 Va Vs V6 Vi 
Vo=[0;0;0]; Vi=[1;0;0]; V2=[1;1;0]; V3=[0;1;0]; Va=[0; 151]; Vs=[0;0;1]; Vs=[1;0;1]; 
V-=[1;1;1]; 


4. SIMULATION 

Table 2 shows in details the DFIG parameters used in this simulation. In this section, a 7.5K W 
DFIG have been simulated by MATLAB/Simulink. In order to observe the regulation behavior of our system 
Figure 1, two main tests have been applied, a) fixed wind profile operation mode and b) aleatory wind profile 
with MPPT control operation mode. 


Table 2. Simulation settings 
Simulated DFIG parameters 


Rated power 7.5kw 
hub height 43,5m 
Moment of inertia 0.3125 kg.m?* 
Gearbox Gain 90 


L, = 0.084H, R, = 0.455 Q , L,=0.081H,R, = 0.62Q,M =0.078H, P = 2, f = 50Hz 


4.1. Fixed wind profile operation mode (test a) 

In this first test, a fix wind speed is applied to the wind turbine blades. It generally consists on 

applying active and reactive power steps that will help us checking the decoupling between the DFIG two 
powers. 
Figure 11 and Figure 12 shows that active and reactive power follows perfectly the reference ones. For 
Figure 13 the sinusoidal waveform of the stator currents is presented. In Figure 14 depicts the variation of the 
mechanical speed versus time. It tracks the wind but with a delay which is due the very fast variation of the 
wind against the slow dynamics of generator. Figure 15 shows the curve of the electromagnetic torque of the 
DFIG. The latter follows the variation of active power. 


Active Power (W)} 
Reactive Power (VAR) 











i L L 
0.5 1 1.5 2 2.5 3 3.5 4 
Time (s) 





Figure 11. Profile of the stator active power Figure 12. Profile of the stator reactive power 
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Figure 15. Time variation of the electromagnetic torque 


4.2. Aleatory wind profile with MPPT control operation mode (test b) 

In the results of the second test, the reference of the reactive power Figure 17 is set zero to ensure a 
unit power factor in the stator side and to optimize the quality of the energy delivered to the grid. While the 
reference of the active Figure 18 power makes it possible to keep the power coefficient of the wind turbine 
optimal whatever the wind speed Figure 16, these results are obtained thanks to the use of the MPPT 
technique. In Figure 19 displays the waveform of stator currents, we notice that the current magnitude varies 
as predicted, and the current frequency is reasonable. In Figure 20 exhibits that the mechanical speed follows 
the desired value computed based on the wind speed. With significant variations around the reference. This is 
because of the fast variations of the wind when we compare it to the generator dynamic. For the (Figure 21) 
shows, the electromagnetic torque of the generator varies according to that of the turbine that is a function of 
the wind speed. 
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Figure 16. Wind speed profile Figure 17. Profile of the stator reactive power 
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Figure 18. Profile of the stator active power Figure 19. Stator current versus time 
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CONCLUSION AND FUTURE WORK 


In this paper, we have presented the modeling of a wind turbine and a doubly fed induction 


generator DFIG. Then, the DPC control strategy has been explained in detail. After that, the simulation 
results demonstrate the effectiveness of the DPC technique. One can notice that the power ripples are 
bounded and the generated power tracks the reference one for both the step the aleatory wind profiles. As a 
perspective of this work, we will compare the performances of DPC classic with those of fractional order 
sliding mode based DPC control by simulation and hardware implementation. 
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